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Abstract 


Tropical glaciers are sensitive indicators of climate change. The loss of the 
volume of tropical glaciers in the Parón basin (Cordillera Blanca, Peru), is 
an example of this at a global level and these changes are expected to 
affect ecosystem services that are unique and irreplaceable. This study 
estimated the economic value of ecosystem services associated with 
water resources impacted by glacial retreat in the Parón basin-Huascarán 
National Park (Cordillera Blanca, Peru) between the years 2009 and 2018. 
The methodology started by mapping glaciers from 2009 to 2018 using 
high-resolution satellite images taken during the dry season. Then, the 
3D surfaces of the glaciers were reconstructed using the GLABTOP tool to 
calculate the volume. With the data obtained, the value of the Tourism 
and Recreation ecosystem service, provision of Water Flow, was 
calculated. The results indicate that the glaciers of the Parón basin have 
been reduced by “1.82% of the glacial area and at a maximum 
deglaciation rate of 0.08 (km?/year). The results of the 3D reconstruction 
show 1.8% more thickness between 40 and 60 m with an ice loss of 1.02 
km3, at a maximum rate of 0.14 km3/year. The present monetary values 
of the services from 2009 to 2018 were estimated at $52 029,34 for the 
tourism and recreation service and $3 213 258.21 for the water flow 
provision service, calculated at 8% discount rate. Therefore, it is 


recommended that investment in programs on the Parón watershed water 
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network, such as increased services, population or industrial use, should 


be increased. 


Keywords: Cryosphere, ecosystem services, economic valuation, 


Cordillera Blanca, climatic change. 


Resumen 


Los glaciares tropicales son indicadores sensibles del cambio climático. La 
pérdida del volumen de los glaciares tropicales en la cuenca Parón 
(Cordillera Blanca, Perú) es una muestra de esto a nivel global y se espera 
que tales cambios afecten a los servicios ecosistémicos que son únicos e 
insustituibles. Este estudio estimó el valor económico de los servicios 
ecosistémicos asociados con el recurso hídrico impactados por el 
retroceso glaciar en la cuenca de Parón-Parque Nacional Huascarán 
(Cordillera Blanca, Perú) entre los años 2009 y 2018. La metodología 
inició cartografiando los glaciares de 2009 a 2018 usando imágenes 
satelitales de alta resolución tomadas en temporada seca. Luego, se 
reconstruyeron las superficies 3D de los mismos usando la herramienta 
GLABTOP permitiendo calcular el volumen. Con los datos obtenidos se 
calculó el valor del servicio ecosistémico del turismo y recreación, y 
provisión de flujo hídrico. Los resultados indican que los glaciares de la 
cuenca Parón se han reducción de “1.82% del área glaciar y en una tasa 
de deglaciación máxima de 0.08 (km?/año). Los resultados de la 
reconstrucción 3D muestran un 1.8% más de espesores entre los 40 y 60 
m, con una pérdida de hielo de 1.02 km?3, a una tasa máxima de 0.14 
km3/año. Los valores monetarios presentes de los servicios de 2009 a 


2018 se estimaron en $52 029.34 para el servicio de turismo y recreación 


2023, Instituto Mexicano de Tecnología del Agua. Tecnología y ciencias del agua, ISSN 2007-2422, 


Open Access bajo la licencia CC BY-NC-SA 4.0 14(6), 76-121. DOI: 10.24850/j-tyca-14-06-03 
(https: //creativecommons.org/licenses/by-nc-sa/4.0/) 


a 1) Check for updates 
OPEN ACCESS 


Tecnología y 


CienciaszAgua 


y de $3 213 258.21 para el servicio de provisión de flujo hídricos, 
calculados a 8% de tasa de descuento. Por lo tanto, se recomienda que 
se debe incrementar la inversión en programas sobre la red hídrica de la 
cuenca de Parón, como mayores servicios, usos poblacionales o 


industriales. 
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Introduction 


Glaciers are freshwater reserves throughout the year; since in tropical 
glaciers ablation is throughout the year (Cogley et al., 2011) providing 
water continuously according to the hydrological cycle (Rumbaur et al., 
2015). The progressive retreat of glaciers in the Andes has been known 
since the twentieth century to the present, and the change in their mass 
balance is a problem that generates great socio-environmental impact 
(Mark €: Seltzer, 2005). According to Urrutia and Vuille (2009) by the year 
2050, all glaciers located below 5 500 masl would disappear, having a 
significant effect on the regulation of the hydrological regime (Jia et al., 


2020). Peru has the largest area of tropical glaciers globally, 
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approximately 71% covering an area of 1200 km2. Most of its population 
resides on the Pacific Ocean coast, one of the driest regions on Earth. Its 
water supply depends largely on glaciers in the Andes Mountains (Kaser 
8 Osmaston, 2002), and to a greater extent in the dry season 
(Schauwecker et al., 2017). The Cordillera Blanca is the most extensive 
glacial mountain range in the tropics (Kaser, Ames, 8 Zamora, 1990)). 
Perhaps for this reason and because of its relative ease of accessibility to 
different glacial areas, it is the most studied area with the greatest 
amount of data available. It is located in the province of Ancash and had 
about 348 km? (Rabatel et al., 2012). Studies in recent years have 
indicated that during 1980 and 2010 the glaciers of the tropical Andes lost 
on average between 0.6 and 1.2 m/year (Rabatel et al., 2013) of their 
mass. But 2004-2014 the Cordillera Blanca presented a glacial retreat of 


31 km?, representing a loss of 6% (Yap, 2015). 


A large number of methodologies have been postulated to 
reconstruct current glacier thicknesses (and indirectly calculate the water 
store), because making direct measurements on the glacier is not possible 
in all mountains (Farinotti et a/., 2017). These numerical modelings can 
be of easy application (Chen £ Ohmura, 1990; Bahr, Meier, € Peckham, 
1997; Lúthi, Walter, Jouvet, 8 Werder, 2006; Radié 8 Hock, 2011; 
Grinsted, 2013) or of higher complexity (Farinotti, Huss, Bauder, Funk, 8 
Truffer, 2009; Linsbauer et a/., 2009; Li, Ng, Li, Qin, 8 Cheng, 2012; 
Linsbauer, Paul, 8 Haeberli, 2012; Farinotti, King, Albrecht, Huss, 4 
Gudmundsson, 2014; Frey et al., 2014; James 8 Carrivick, 2016). The 
calculation of volume in the Peruvian Andes has not been widely 


estimated. Currently there are data for point glaciers in the Cordillera 
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Blanca, as well as for areas south of the Andes. The most commonly used 
methodology so far is GlaBTop (Linsbauer et a/., 2009; Linsbauer et al., 
2012; Iparraguirre et al., 2020) for the estimation of future potential risk 
lakes but which also allows the calculation of the water store (Iparraguirre 
et al., 2020). 


Glacier retreat and anomalies in precipitation patterns, pose a 
serious threat in water availability (Bury et a/., 2011) and high variability 
that will significantly affect ecosystem services (Heikkinen, 2017) being a 
potential source of social and political tensions that can turn into 
instability or even serious conflicts (Hijioka et a/., 2014). These conflicts 
could result in restricted access to certain services such as the use of 
potable water for drinking and irrigation, heavy dependence on 


agriculture, and limited alternative livelihood opportunities (Das, 2009). 


Evaluations of these services can be performed through economic 
valuation, using market value-based methods, revealed preference 
methods, and stated preference methods (MINAM, 2016). Currently, 
there are few studies of its application in China (Yuan 8 Wang, 2018; 
Zhang et al., 2019), Spain (Grima 8 Campos, 2020), Chile (Segovia, 
2014) and Peru (Boyano, 2016). Therefore, identifying and assessing the 
impacts generated by glacial retreat on ecosystem services will allow 
decision making for conservation, prevention measures, control and 


development of targeted projects (Boyano, 2016). 


The Parón-Huascarán National Park (Cordillera Blanca) watershed 
has its own intrinsic characteristics such as the paramos (lower part of 
the mountain ranges) and the Cryosphere (upper part of the mountain 


range) that provide ecosystem services for the populations in the 
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watershed (Boyano, 2016) such as support, provisioning and regulation 
(Liekens et al. 2013). Glaciers in the Parón basin are highly vulnerable to 
climate change, which have induced significant changes in volume loss, 
surface area and permafrost loss (Medina 8: Mejía, 2014; Navarro et al., 
2017) as already discussed. However, this vulnerability directly affects 
ecosystem services related to water resources and there is currently no 


quantification of their value. 


Therefore, this research attempts to fill this gap by (1) examining 
the changes in glacial retreat as surface area and volume lost between 
the years 2009 and 2018; (2) approximating to obtain a baseline value of 
the ecosystem services of: tourism and recreation and continuous water 


flow provision, using the real markets method for this purpose. 


Methodology 


Study area 


The Cordillera Blanca is located in northern Peru and is the most extensive 
glaciated mountain range in the tropics (Kaser et al., 1990). For all 
analyses in this study, we refer to the Parón watershed, one of the ten 
glaciated watersheds within and adjacent to the Callejón de Huaylas 
(Figure 1). 
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Figure 1. Location of the glacial basins of the Cordillera Blanca Source: 
Mark and Seltzer (2003). 


The Parón basin is located north of the Cordillera Blanca +90 km 


east of the Pacific Ocean (taking Isla Blanca del Callao as a reference). It 
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is located in the district of Caraz, province of Huaylas, in the department 


of Ancash, “57 km from the city of Huaraz, extending between 8%S/77W., 


It drains its waters into the Santa River basin and one of the first 


inventories estimated that in 1962, 55% of its area was glaciated (Table 


1). 


Table 1. Physical characteristics of the hydrographic basins of the 


Cordillera Blanca. 


Glaciated , Average Min. Max. 
Basin area A Ñ E 
Watershed area (1962) (km2) elevation elevation elevation 
m 
(%) (aslm) (aslm) (aslm) 
Parón 55 41 4920 4152 5965 
Llanganuco 41 89 4832 3831 6670 
Marcara 
25 259 4454 2887 6196 
(Chancos) 
Cedros 22 114 4544 1975 6128 
Colcas 19 236 4313 2022 6178 
Quilcay 18 240 4522 3131 6195 
Pachacoto 12 206 4610 3713 5574 
Olleros 11 174 4437 3461 5673 
Santa (La Balsa) 9 4784 4056 1858 6733 
Quitarasca 8 384 4237 1587 5921 
Querococha 6 63 4524 3992 5291 


Source: Taken from Mark and Seltzer (2003). 
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Tourism in glacial areas is also an important source of income and 
employment for local communities adjacent to tourist sites (Haimayer, 
1989). The Parón basin is an important tourism and recreation site within 


Huascarán National Park. 


Extraction of information for the change of the glacier 
area 


The delimitation of the glacier surface, as well as the 3D reconstruction of 
its surfaces was carried out using ArcGIS 10.5 software. The materials 
used as inputs were: satellite images obtained from Google Earth in dry 
seasons corresponding to each year of study (from 2009 to 2018); a 
Digital Elevation Model (or DEM) with 30 m resolution corresponding to 
the years of analysis; and, finally, the contour lines (equidistance=50 m) 
of the study area. All these inputs were provided by the Instituto 
Geológico, Minero y Metalúrgico (INGEMMET). 


Delimitation of current glaciers and paleoglaciers 


The delimitation of the glacier extent was based on a Google Earth 


satellite image of 30 meters per pixel resolution. 


The delimitation was performed on the satellite image with the 


ArcMap mapping tool of ArcGIS. Each glacier device was manually 
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digitized from the cirque to the glacier tongue at a mapping scale of 
1:250. In addition, other inputs were used such as the Hillshade of the 
satellite image (30m/px), 1962 aerial photographs of the study area and 
supporting tools such as ArcScene and historical images from Google 
Earth. 


To avoid an overestimation of the glacier area due to a possible 


snow mapping, a low-water season satellite image was chosen. 


Reconstruction of the 3D geometry of glacial and 
paleoglacial surfaces 


Knowing the 3D geometry of current and past glaciers is of great 
importance because it allows us to quantify the mountain's water 
reserves, its deglaciation rate and therefore its evolution over a given 


period of time. 


The 3D reconstructions of the Peruvian Andes have been carried out 
using Ground Penetrating Radar (GPR) and because Andean glaciers are 
very difficult to access, GPR has been used in specific glacier mountains: 
Artesonraju (Cordillera Blanca, Ancash), Coropuna (Cordillera Occidental, 
Arequipa) and Quelccaya (Cordillera Oriental, Cusco-Puno). Therefore, it 
is necessary to complement the GPR work with indirect measurements 
using numerical models that reconstruct the glacier surface as close to 


reality as possible. 
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Currently, the most widely used  methodologies for 3D 
reconstruction of current and past glaciers are GlaBTop and GlaRe, 
respectively. 


GlaBTop 


This methodology allows estimating ice thicknesses along flow lines, 
which represent the path of the glaciers (Figure 2). Subsequently, the 
thicknesses were interpolated (Figure 3) covering the entire glacier area, 
obtaining a raster (3D representation). The sum of all the thicknesses 
multiplied by the cell size of the raster made it possible to obtain the water 
volume represented by the mountain. On the other hand, these 
thicknesses were subtracted from the current Digital Elevation Model 
(original DEM minus thickness DEM), to obtain the basal topography 
(Linsbauer et a/., 2009; Paul € Linsbauer 2012). 
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Figure 2. Mapping of o mes in some glaciers of the Parón basin on 
the Hillshade of the digital elevation model (DEM). 


Figure 3. Interpolation of thicknesses using y the" 'topo to raster" tool of 
the GlaRe tool. 
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GlaRe 


GIS tool (Pellitero et a/., 2016), which semi-automatically reconstructs 
the three-dimensional geometry of paleoglaciers (in this article all 
scenarios for years prior to 2018), based on the basal topography 
calculated with the GLABTOP methodology and new flow lines of its own 
for each scenario. The tool requires three input elements: On the one 
hand, the flow lines converted into centroids (where the generated 
thicknesses will be stored), the DEM of the basal topography, and finally, 
since the model does not reflect the stress generated by the lateral walls 
of the valley, shape factors were required based on evidence of glacier 
thickness observed by satellite images, lateral moraines and valley 
geometry following the equation proposed by Nye (1952), adjusting the 
thicknesses generated in the paleoglacial front. 


Method for estimating the economic value of ecosystem 
services 


To measure the Ecosystem Service Value of Tourism and Recreation, the 
Market Price Method was used, as it is the best known and allows 
estimating direct use values (MINAM, 2015). It is proposed that at least 
the value of the ecosystem service is equal to the value incurred by 
visiting them, it has been established that for each entry fee to Huascarán 


Park,% corresponds to the ice present inside the park, /.e., that amount 
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would correspond to the monetary value of the ecosystem service 
(Segovia, 2014) of tourism provided by the entire surface area of 
Huascarán National Park (Cordillera Blanca) and the equivalent value of 
the glaciers present in the park. The entrance fee in the last 10 years has 
been between S./15.00 (<> 3.87 dollars) for national tourists and S./ 


30.00 (<> 7.74 dollars) for international tourists: 


As = VPse (1) 
As =Ag (2) 
Where: 


As: Surface area of Huascarán Park. 
VPs: Present value of the tourism and recreation service. 


Ag: Glacial area of lake Parón as part of Huascarán National Park. 


For the monetary valuation of the glaciers present in the study area, 
it was established that for each entry value to Huascarán National Park, 
there is a certain amount of surface area of the glaciers present in the 
Parón basin, ¡.e., this amount of glacier would correspond to the monetary 
value of the ecosystem service of tourism and recreation that has been 
lost due to glacier retreat. Therefore, the value of the deglaciation rate 


due to loss of ecosystem service would be equal to the% loss of surface 
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area in the 2009-2018 time periods. 
VTD = % Agp x VPse (3) 
Where: 


VTD: Value of the rate of deglaciation due to loss of ecosystem service. 
VPse: Present value of Tourism and Recreation Service. 


Agp: Glacial area retreated in the Parón watershed from 2009 to 2018. 


The calculation of the monetary value of the ecosystem service of 
"provisioning of continuous water flow" was estimated based on the 
average annual flow expressed in cubic meters per second (m/s) 
provided by the glaciers of the Parón-Huascarán National Park (Cordillera 
Blanca) basin and the price of water rights markets are according to the 
Supreme Decree that approves values of economic retribution to be paid 


for surface water use (El Peruano, 2018): 
VSEryc = Pr x QMac (4) 
Where: 


VSErHc = Continuous water flow ecosystem service value. 
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Pr = Non-consumptive DAA price. 


QMacs = Mean annual glacial inflow (m?/s). 


Estimation of the economic value of glacial retreat 


The monetary value was estimated based on real prices of the ecosystem 
services of "Tourism and Recreation", "Water Storage" and "Continuous 
Water Flow" of glacial retreat in the Parón Basin. For this, the annual 
values (annual economic benefits) were calculated for being within a 
protected site, the ecosystem services valued by the glacier loss will 
maintain such services for a defined period of time. In addition, to obtain 
the present value (PV) of the valued ecosystem services, a social discount 
rate of 8% was used (MINJUS, 2017), which is the one considered for 


social projects by the Peruvian Ministry of Justice: 


Fi=Pix(1+n" (5) 
Where: 
Fi = Initial value of the service flow ¡ (annual economic benefit). 


P¡ = present value of the flows, the stock value of the service. 
T = Time. 


r = Constant discount rate for the whole period (8%). 
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Results and discussion 


Variation of glacier characteristics 


Glacier delimitation 


The behavior of glaciers in the Parón Basin during this 10-year period 
must be inferred from the retreat of glacier surface cover suffered by the 
Cordillera Blanca that has accelerated since 1987 (Burns € Nolin, 2014). 
And considering that the observations of Georges (2004), already 
estimated a 1% increase in the rate of glacier loss per year from 1987 to 
1996 in the Cordillera Blanca. 


The calculation of the extension for each period analyzed in this 
work is presented in Figure 4. 
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Evolution of the glacial area in the Parón-Huaraz basin 


Area km2 


Glacier area estimation 


Figure 4. Comparison of glacier area decrease between 2009-2018 in 


the Parón watershed. 


All these results have only considered discovered ¡ce masses where 
the glacial extent of the Parón Basin for 2009 covered an area of 18.82 
km?, while for 2018 it was reduced to 18.48 km?. This meant a reduction 
of “1.82% of the glacier area and can be translated into a deglaciation 
rate (km?/year) minimum from 0.02 and maximum of 0.08 from 2009 to 
2018; meaning that the deglaciation rate in the last 10 years remains 
constant (Figure 5). 
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Figure 5. Satellite model of glacier cover retreat between 2009-2018. 
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Thus, our estimates of glacial change are aligned with important 
works (Georges (2004), Burns € Nolin, 2014) that agree on an 
acceleration in the loss of glacier surface area in the Cordillera Blanca 
since 1970, where it went from 723 to 482 km?, 


3D reconstruction of the geometry of glaciers and 
paleoglaciers 


Deducing how much the glacier area has been reduced in terms of area, 
in a way, can be misleading. Because when viewed in a 2D plane, we 
cannot perceive how the volume of water in solid state has changed, since 
the thicknesses that would be disappearing are unknown. For this reason, 
this research applies novel numerical modeling that can give an indication 
of how glacier masses are evolving. It is worth mentioning that the 
glaciers delimited, the greatest range of thicknesses oscillate between 20 
and 40 m, but the great difference lies in the prevalence of ice masses of 
greater thickness (Figure 6). That is, the difference in glacier volume in 
2009 was maintained with greater power; whereas, by 2018 they have 
been affected. 
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Figure 6. Satellite 3D reconstruction model of glacier volume loss 
between 2009-2018. 
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Thickness distributions reveal that on the 2018 glaciers (maximum 
thickness = 180 m on snow-capped Artesonraju) thicknesses between 20 
and 40 m predominated. These thicknesses cover 68% of the glacier area 
representing “11 km?. The greatest thicknesses are found specifically in 
the Artesonraju (maximum thicknesses of 180 m and greater than 100 
m). This is explained by the fact that this mountain is considered as an 
ice field; that is, where the slope is low. Glaciers in 2009, although it is 
true, anyway, predominantly thicknesses between 20 and 40 m (66% <> 
12 km), there were still thicknesses between 40 and 60 m (at least 12% 


more compared to 2018). 


For example, Colonia et al. (2017) estimated using this 3D 
reconstruction methodology the possible lake formations throughout the 
Peruvian Andes. This work is an excellent indicator of the areas with 
greater thickness at the regional level, but at the local level it loses 
precision. One of the drawbacks is that the input DEM has a 90 m 
resolution. Torres et al. (2014) calculated the total volume of the Peruvian 
Andes based on two approaches: area-volume and slope-thickness being 
35 and 34.39 km, respectively. These results have also been based on 
low resolution DEMs (60 m) that are useful for large-scale estimates, but 
do not allow obtaining a detailed basal topography necessary for 
paleoglacial reconstruction. But at the local level, compared to those 
estimated in this work, the results are close. Basically, these small 


differences are linked to the processing and input data (Figure 7). 
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Figure 7. Glacier volume reduction in the Parón basin 2009-2018. 


From the glaciers of 2009, we have 1.8% more thicknesses between 
40 and 60 m compared to the glaciers in 2018 with 1.02 km? lost, with a 
maximum rate of 0.14 km?*/year. This difference between these scenarios, 
show us that, despite having a low deglaciation rate in reference to the 
estimated area, we can deduce that the best indicator is the 3D 


reconstruction because it allows us to observe what potential ¡ce volume 
has been lost in that area. 
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Ecosystem service value analysis associated with 
glacial retreat 


Tourism and recreation ecosystem service valuation 


Tourism in the Cordillera Blanca, the highest mountain range in the 
Peruvian Andes, has developed greatly in the last 30 years due to the 
number of domestic and foreign tourists visiting the area (Grótzbach, 
2003). Huascarán National Park is the largest park in the Cordillera 
Blanca, promoting many local activities and influencing environmental, 
social, economic and tourism aspects of the region (Maguiña, Angulo, 
Gonzales, 8: López, 2020). 


The total number of visitors to Huascarán National Park in the last 
10 years (2009-2018) was 990 306 for domestic tourism and 320 065 
foreign tourists, showing an annual growth trend (Figure 8a). There are 
more than 100 known mountain glacier resorts in the world that operate 
as special natural attractions (Wang, Zhao, € He, 2012); among the most 
important ones the Huascarán National Park that in the years after 2014, 
glacier tourism in the different basins has been developed and operated 


successfully, obtaining high economic benefits (Izaguirre, 2021). 
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Figure 8. Variation in tourism and economic value of glacial retreat in 
the Parón Basin (2009-2018). Left (a): Number of visitors per year. 


Right (b): Variation of monetary income as a function of glacial retreat. 


With the rapid growth of the world economy, accessibility to 
transportation and the increasing demand of visitors tourism has 
improved and increased considerably, positioning itself with an average 
of 163 796.375 visitors per year, being an average income compared to 
the world trend of other high mountain glacier tourist sites as an example 
the national park of "Glacier Bay" (Alaska) with 400 000 visitors per year 
(Pirhalla, Gende, 8: Mólders, 2014); Athabasca glaciers (Canada) with 600 
000 visitors per year; Glaciares national park (Argentina) with 167 000 
visitors per year; "Jostedal Glacier" national park (Norway) with 40 000 


visitors per year (Welling, Árnason, 8 Olafsdottír, 2015), and "Franz 
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Jozef" Glacier (New Zealand) with 346 000 visitors per year (Purdie, 
2013). 


Figure 8b shows that economic value is being lost whenever there 
is a loss of glacial area in Parón, it is recognized that glacial retreat 
significantly impacts tourism and recreation (Vuille et a/., 2018) reducing 
the number of visitors in tourism-dependent cities (Purdie, Gomez, € 
Espiner, 2015), progressively affecting the economic value of the Parón 


watershed that it can provide as an ecosystem service. 


The valuation method of personal preferences based on real prices 
used to calculate the monetary value of the ecosystem service of tourism 
and recreation showed that the total equivalent to the surface area lost 
due to glacial retreat between 2009 and 2018 in the interior of Huascarán 
National Park has a total economic value of s. /133 788.19 (equivalent to 
$33 122.46 dollars at an exchange rate of 4.04), which, updated to 
present value at a rate of 8%, corresponds to a value of S/210 156.82 
(equivalent to $52 029,34) (Table 2). This shows that there is an annual 
increase in monetary value lost to glacial retreat. Studies such as Zhang, 
Zhuoran and Seenprachawong (2015) estimated that annual monetary 
valuation of glaciers in Mt. Yulong, Yunnan, China for tourists was $500 
million, while the studies conducted to the El Morado glacier in Chile 
(Segovia, 2014) shows monetary values of $25 554 annually, with values 
closer to those obtained in this study, indicating that the value of glaciers 


for the tourism sector is substantial. 


Table 2. Annual flows of the monetary amount of the ecosystem service 


of Tourism and Recreation in the Parón watershed. 
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Total annual 
Total annual e 
deglaciation value 
monetary value Surface 
N* | Year | National | Foreigners , (S./ /m?/year) 
to Parque area loss 
£ Monetary Current 
Huascaran 
value value 
1 |2009| 66278 31071 S/1,926,300.00 | 18819080.00 | S/10,662.12 | S/11,589.26 
2 |2010| 79852 32067 S/2,159,790.00 | 18755666.21 | S/11,914.21 | S/14,076.33 
3 |2012 | 93635 33185 S/2,400,075.00 | 18704096.37 | S/13,203.30 | S/16,955.83 
4 |2014 | 103584 33950 S/2,572,260.00 | 18686107.44 | S/14,136.92 | S/19,733.47 
5 |2015| 112818 35758 S/2,765,010.00 | 18604096.37 | S/15,129.56 | S/22,955.53 
6 |2016 | 139063 38799 S/3,249,915.00 | 18538561.70 | S/17,720.22 | S/29,224.17 
7 |2017| 200189 48971 S/4,471,965.00 | 18500997.58 | S/24,334.06 | S/43,621.42 
8 |2018 | 194887 66264 S/4,911,225.00 | 18475738.22 | S/26,687.80 | S/52,000.81 


Valuation of the Ecosystem Service of Water Flow 
Provisioning 


Deglaciation as part of the natural dynamics of the Cryosphere a 


provisioning function to meet a social-ecological demand and may pose 


new challenges for societies (Zhang et a/., 2021b). The use of water flow 


provision due to deglaciation 


is not only reflected 


in an annual 


contribution, but also in a seasonal one, which is currently undergoing 


alterations due to climate change. 


The water flow provision service is subject to an increase in 


temperatures due to climate change, which, generates a surplus of 
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snowmelt increasing the water flow in the basin and a decrease of 
permafrost in the ¡ce sheet (Sotillo, 2017). Permafrost is scarce in the 
Cordillera Blanca; therefore, melting due to temperature increase causes 
an increase in the water flow provision service in the different basins 
(Iparraguirre et al., 2020), thus giving rise to unusual water flow 


increases. 


Table 3 shows that the economic value of the ecosystem service 
provision of water flow due to the natural dynamics of the permafrost 
added to the deglaciation due to climate change in the Parón watershed 
in the last 10 years was S./7 762 276.36 with a present value of S./12 
981 563.18 (equivalent to $3 213 258.21 with an exchange rate of 4.04) 


at an 8% discount rate. 


Table 3. Annual flows of monetary value foregone for the ecosystem 


service of hydrological flow provisioning by type of use in the Parón 


watershed. 
Economic Total annual deglaciation 
Flows discharged - 
É > Type of surface retributions to be value (S./ /m*/año) 
N2 | Year by glaciers E . 
me vean) water use right paid for water use Monetary Current 
3 
mo value value 
Poblacional S/0.0050 S/22,826.50 | S/24,811.41 
Industrial* S/0.0792 S/361,571.75 | S/393,012.77 
id cid Minero 3/0.1017 5/464,290.99 | 5/504,664.12 
Other uses ** S5/0.0330 S/150,654.89 | S/163,755.32 
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Flows discharged 


Type of surface 


Economic 


retributions to be 


Total annual deglaciation 
value (S./ /m*/año) 


N2 | Year by glaciers . . 
(me/yeao) water use right paid for water use Monetary Current 
3 
(Uno value value 

Poblacional S/0.0050 S/22,784.49 | S/26,919.29 

Industrial* S/0.0792 S/360,906.30 | S/426,401.59 
2 2010 4556897.77 

Minero S/0.1017 S/463,436.50 | S/547,538.40 

Other uses ** S/0.0330 S/150,377.63 | S/177,667.33 

Poblacional S/0.0050 S/22,779.44 | S/31,797.41 

Industrial* S/0.0792 S/360,826.39 | S/503,671.05 
4 2012 4555888.81 

Minero Ss/0.1017 S/463,333.89 | S/646,759.42 

Other uses ** S/0.0330 S/150,344.33 | S/209,862.94 

Poblacional S/0.0050 S/22,578.80 | S/37,236.94 

Industrial* S/0.0792 S/357,648.26 | S/589,833.12 
6 2014 4515760.89 

Minero Ss/0.1017 S/459,252.88 | S/757,399.35 

Other uses ** S/0.0330 S/149,020.11 | S/245,763.80 

Poblacional S/0.0050 S/22,379.93 | S/40,118.43 

Industrial* S/0.0792 S/354,498.12 | S/635,475.99 
7 2015 4475986.41 

Minero Ss/0.1017 S/455,207.82 | S/816,008.95 

Other uses ** S/0.0330 S/147,707.55 | S/264,781.66 

Poblacional S/0.0050 S5/21,665.14 | S/42,214.22 

Industrial* S/0.0792 S/343,175.75 | S/668,673.28 
8 2016 4333027.11 

Minero s/0.1017 S/440,668.86 | S/858,637.28 

Other uses ** S/0.0330 S/142,989.89 | S/278,613.87 
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Economic Total annual deglaciation 
Flows discharged A Ea 
- Type of surface retributions to be value (S./ /m*/año) 
N2 | Year by glaciers . . 
(me/yeao) water use right paid for water use Monetary Current 
3 
(Uno value value 
Poblacional S/0.0050 S/21,335.41 S5/45,186.69 
Industrial* S/0.0792 S/337,952.86 | S/715,757.12 
9 2017 4267081.52 
Minero S5/0.1017 S/433,962.19 | S/919,097.21 
Other uses ** S/0.0330 S/140,813.69 | S/298,232.13 
Poblacional S/0.0050 S/20,952.17 | S/48,233.72 
Industrial* S/0.0792 S/331,882.36 | S/764,022.07 
10 2018 4190433.90 
Minero s/0.1017 S/426,167.13 | S/981,073.79 
Other uses ** S/0.0330 S/138,284.32 | S/318,342.53 


Studies such as Jeong, Sushama and Naveed-Khaliq (2017) show 
that in the Northern Hemisphere the average snow cover mass decreases 
at a rate of 10.2 Mt per year, or decreases by 0.3% per year due to a 
global warming effect. So also, projecting to 2030 the effects of climate 
change, the value of snowpack losses in western China would be reduced 
by CNY2.9 billion (Wu et a/., 2021). This research suggests that climate 
change is affecting all glaciers worldwide and will have a direct effect on 


the economic values of ecosystem services. 


Direct effects such as those shown in this study are also present in 
glaciers in the Quilian regions (China) where they show that the service 


value of glacier meltwater provision to river runoff was 8 072 million yuan 
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equivalent to 33% ($1 277 million USD), the service value of freshwater 
supply was 844 million yuan ($133 million USD) (Sun et a/., 2020). 


Some studies such as in Latin America as Segovia (2014) has as a 
result that the present value of the ecosystem service of continuous water 
flow from the glaciers of El Morado Natural Monument, Chile, is clp$2 327 
683 926 ($15.53 million dollars) and the annual value of this service with 
a discount rate of 6% reaches clp$139 661 036 ($175 thousand dollars). 
This study is the closest to the reality present in the Parón basin, the 
economic values are not as far away as those presented by Wu et al. 
(2021) and Sun et al. (2020). 


Regarding the service of water flow provision, it is necessary to 
emphasize that the cost of water is based on market prices for its type of 
use. In general, the cost of water use in Peru is reflected in terms of 
mining, urban, rural and industrial use. As shown in Figure 9, the present 
value of the ecosystem service of water flow provision in the industrial or 
mining sector has the highest values as opposed to the other types of 
use. Similarly, the study by Zhang et al. (2021b) reports that the average 
annual economic value of water uses for agriculture, services, urban use 
and industry estimated at 1.7, 4.4, 9.9 and 15.7 million yuan with the 


industrial sector having the highest consumption. 
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Figure 9. Variability of the economic value of the ecosystem service of 
water flow provisioning according to the type of use in the Parón 


watershed. 
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Conclusions 


This research presents the glacier loss in the Parón Basin (Cordillera 
Blanca), the retreat of which directly influences the ecosystem services 
associated with water resources. The Parón watershed is characterized by 
the tourism/recreation service and the provisioning of water flows 
presenting a high economic valuation that has grown in recent years 
(2009-2018). 


The Parón basin in 2009 covered an area of 18.82 km?, while by 
2018 it was reduced to 18.48 km?. This meant a reduction of “1.82% of 
the glacier area and can be translated into a deglaciation rate (km?/year) 
minimum from 0.02 and maximum of 0.08 from 2009 to 2018. Similarly, 
from the 2009 glaciers, we have 1.8% more thicknesses between 40 and 
60 m compared to the glaciers in 2018 with 1.02 km3 lost, with a 
maximum rate of 0.14 km3/year. These results evidence an effect of 


climate change on the Parón watershed. 


The economic value of the ecosystem service of tourism and 
recreation was calculated using the method of valuation of personal 
preferences based on real prices. The results show that the monetary 
amount lost due to the retreat of the glacier surface area between 2009 
and 2018 in the interior of Huascarán National Park is s./133 788.19 
(equivalent to $33 122.46 dollars at an exchange rate of 4.04), which, 
updated to present value at a rate of 8%, corresponds to a value of S/210 
156.82 (equivalent to $52 029.34). 
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The corresponding economic value of the water flow provisioning 
services in the different water uses (Population, Industrial, Mining, Other 
uses) generated by the meltwater from the Parón watershed between 
2009 and 2018 was S./7 762 276.36 with a present value of S./12 981 
563.18 (equivalent to $3 213 258.21 with an exchange rate of 4.04) at a 
rate of 8%. 


Although it was complex to capture the values of the services of 
glacial retreat in the basin in economic terms, it is shown that the service 
of water flow provision has a higher monetary value than recreation and 
tourism, suggesting that climate change associated with deglaciation is 
generating a loss of water value that is not being taken advantage of. It 
is recommended that current investment in programs on the water 
network of the Parón watershed should be increased, such as greater 


services, population or industrial use. 
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